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Through reformulating the cold, self-gravitating fermion gas as a Bose condensate by identifying 
their mutual Thomas-Fermi limits, the dissipationless formation of a heavy neutrino star in gravi- 
tational collapse is numerically demonstrated. Such stars offer an alternative to supermassive black 
holes for the compact dark objects at the centers of galaxies. 
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Supermassive neutrino stars, in which gravity is bal- 
anced by the degeneracy pressure of cold fermions, have 
been a subject of speculation for more than three decades 
0. Originally, these objects were proposed as models 
for dark matter in galactic halos and clusters of galaxies, 
with neutrino masses in the ~ eV range. More recently, 
however, degenerate superstars composed of weakly in- 
teracting fermions in the ~ 10 keV range have been sug- 
gested as an alternative to the supermassive black holes 
that are purported to exist at the centers of galaxies |§-D . 
In fact it has been shown Q that such degenerate fermion 
stars could explain the whole range of supermassive com- 
pact dark objects which have been observed so far, with 
masses ranging from 10^ to 3 x IQ'^Mq, merely assum- 
ing that a weakly interacting quasi-stable fermion of mass 
TO/ ~ 15 keV exists in nature. 

As an example, the most massive and violent compact 
dark object ever observed is located at the center of M87 
with a mass M ~ 3.2 x lO^M© 0. Interpreting this 
as a relativistic fermion star at the Oppenheimer-Volkoff 
1^ limit yields the fermion mass to/ ~ 15 keV and a 
radius R = 4.45 i?s ~ 1-5 light-days ||,||, where Rs is 
Schwarzschild radius. In this case there is little difference 
between the fermion star and black hole scenarios be- 
cause the last stable orbit around a Schwarzschild black 
hole is at 3i?s anyway. 

Extrapolating this down to the compact dark object 
at the center of our galaxy with mass M ~ 2.6 x 
lO^M©, which is at the lower limit of the mass range, and 
i? ^ 20 light-days, the same fermion mass gives R ~ 
10^ i?s- Consequent upon the shallow potential inside 
this fermion star, the spectrum of radiation emitted by 
accreting baryonic matter is cut off for frequencies larger 
than 1013 Hz |,|, as is observed in the spectrum of the 
strong radio source Sgr A* at the galactic center [0. 
This fermion star is also consistent |^ with the observed 
motion of stars within a projected distance of 10 to 30 
light-days of Sgr A* §. 

Of course, it is well-known that 15 keV lies squarely in 
the cosmologically forbidden mass range for stable active 
neutrinos v pd] ]. Sterile neutrinos are another matter: as 
shown by Shi and Fuller 1 12|, in the presence of an initial 
lepton asymmetry of ~ 10^ , a sterile neutrino of mass 



TOs ~ 10 keV is resonantly produced with near closure 
density, $1 = 1. Moreover, the resulting energy spectrum 
is not thermal but rather cut off so as to approximate a 
cold degenerate Fermi gas. This model is constrained by 
astrophysical bounds on |p^ , however the al- 

lowed parameter space includes TOs ~ 15 keV contribut- 
ing ~ 0.3 as favoured by the BOOMERANG data 

The statics of degenerate fermion stars is well under- 
stood, being the Oppenheimer-Volkoff equation in the 
relativistic case [|| or the Lane-Emden equation with 
polytropic index n = 3/2 in the nonrelativistic limit 
Alternatively, because R 3> l/rnf, one may understand 
these as the Thomas-Fermi theory applied to self- gravi- 
tating systems. The extension of the Thomas-Fermi the- 

has been used to show 



ory to finite temperature |16 
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that at a certain critical temperature weakly interacting 
massive fermionic matter undergoes a first-order gravita- 
tional phase transition from a diffuse to a clustered state, 
i.e. a nearly degenerate fermion star. Such studies do 
not, however, bear on the crucial dynamical question of 
whether the fermion star can form through gravitational 
collapse of density fluctuations in an orthodox cosmolog- 
ical setting. Indeed, since collisional damping is negligi- 
ble, one would expect that only a virialized cloud results 

iV-body simulations of the coUisionless Boltzmann or 
Vlasov equation evidence a rather different picture: the 
collapse is followed by a series of bounces with matter 
expelled at each bounce, leaving behind a condensed ob- 
ject |l8|. By Liouville's theorem the Vlasov equation 
describes an incompressible fluid in phase-space so that 
it respects a form of the exclusion principle. Hence, 
these A^-body simulations are effectively fermion simu- 
lations. What transpires is that gravity, being attrac- 
tive, self- organises the phase-space fluid into a high- 
density/momentum core at the expense of other low- 
density/momentum regions as seen in the evolution of 
the spherical Vlasov equation [|l9| . 

Much the same behaviour is observed in the forma- 
tion of mini-boson stars through so-called gravitational 
cooling ||2^. Such a mini-boson star is stable by bal- 
ancing uncertainty and gravitational pressure. A simi- 
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lar mechanism works in the presence of a quartic self- 
interaction which dominates over uncertainty pres- 
sure resulting in an equilibrium radius 1 /rrib where 
nib is the boson mass Hence we have a universal 

description of the physics underlying the formation pro- 
cess: once the collapse proceeds far enough (uncertainty, 
interaction or degeneracy) pressure results in a bounce, 
the outgoing shock wave carrying away the binding en- 
ergy. The virial argument above is circumvented because 
the ejected matter invalidates its assumption that there 
is no flow through the boundary. 

In this letter we verify the above picture for the for- 
mation of the fermion star from a cold nonequilibrium 
configuration. The dynamical Thomas-Fermi theory was 
given long ago by Bloch for the electron gas and 
amounts to Euler's equations for irrotational flow to- 
gether with an equation of state P = P{p)- The prob- 
lem is that, transcribed to the self-gravitating fermion 
gas, there is the Jeans instability, signalled by an imagi- 
nary plasma frequency, and thus short-wavelength shocks 
must be regulated. The usual remedy is to introduce 
some small numerical viscosity; however it seems impru- 
dent to draw conclusions based on introducing dissipa- 
tion into what is fundamentally a dissipationless process. 
Here we take another, literally conservative approach. 

In the Newtonian limit a self-interacting boson star is 
governed by the Gross-Pitaevskii-like equations 
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where for convenience we have absorbed the boson mass 
TOfc in the field. Using the ansatz ij) = ^/p exp ( - irubO), 
we arrive at 
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The Thomas-Fermi limit is governed by nib 3> | V/o| /p. 
Thus neglecting the last term in (2b), Bloch's equations 
are recovered, with being the velocity potential and 
V{p) given by 
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For a general polytropic equation of state 
P{p) = Xpi+i/" 

eq.(||) yields 



(3) 



(4) 



V{p) = (n + 1) if mb pi/" . 
Using the potential (|^) and introducing 
A _[{n + \)Kf''^ 
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as the length and mass scales, respectively, the substitu- 
tion 



yields the dimensionless equations 



A 'dt 
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The validity of the Newtonian approximation in the static 
case requires 

M/A/, = (47r)~i / (fr < 1 , n < 3 . (9) 



For weakly interacting degenerate fermions, the poly- 
tropic index is n ~ 3/2, and 
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0.2325 f— V W-lyr, 
(10a) 



= 1.185 
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where 5/ is the spin degeneracy and mpi is the Planck 
mass. By construction a large but finite A allows us to 
simulate the fermionic problem as a bosonic one through 
their mutual Thomas-Fermi limits, while providing an ex- 
plicitly energy conserving way of controlling the shocks. 
The basic regulating mechanism is the "kinetic" part of 
( ^ ) which penalises gradients of order A. Of course, A 
must be sufficiently large that this term does not change 
the static scaling relationship 



M = const 



(11) 



arising from the polytropic equation of state. Our crite- 
rion is that the ratio of "kinetic" and "pressure" contri- 
butions to the static energy functional should be small, 
in particular for a Gaussian \1/ = a exp [— (r//3)^] 
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FIG. 1. Contour plot for the evolution of jr^'p from the 
initial configuration in text. 



weakly interacting fermionic matter will form supermas- 
sive fermion stars through gravitational collapse accom- 
panied by ejection. For a fermion mass of rn/ ~ 15 keV 
such a superstar is consistent with observations of the 
compact dark object at the center of our galaxy. A sim- 
ilar demonstration for formation near the Oppenheimer- 
Volkoff limit, and the question of cosmology with degen- 
erate dark matter requires a general relativistic exten- 
sion which is under development and will be reported 
elsewhere. 
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FIG. 2. Ir*!"! versus r on the final time slice in Figure 1. 

In Fig. 0we display the evolution of \r^!\'^ for spherical 
collapse of a mass M = 0.008 Af,, initially in the form 
of a Gaussian with /? = 100 = A. The expected features 
of bounce and ejection leaving a condensed core are ev- 
ident. Here we have implemented a velocity dependent 
imaginary part to the potential in the outer layers of the 
cavity to remove the ejected fermion matter before it can 
be artificially reflected by the boundary. The core size 
r/i?* = 26 and mass m/M* = 0.0057 is commensurate 
with a fermion star, however it is far from smooth as 
evidenced by the plot of \r'^\'^ on the final time slice in 
Fig. ||. This feature may, however, be attributed to the 
relatively short duration of the simulation. 

In summary, using a bosonic representation of the 
dynamical Thomas-Fermi theory for a self-gravitating 
gas, we have shown that nonrelativistic, degenerate and 
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